Statistics of Heat Transfer in Mesoscopic Circuits 
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A method to calculate the statistics of energy exchange between quantum systems is presented. 
The generating function of this statistics is expressed through a Keldysh path integral. The method 
is first applied to the problem of heat dissipation from a biased mesoscopic conductor into the ad- 
jacent reservoirs. We then consider energy dissipation in an electrical circuit around a mesoscopic 
conductor. We derive the conditions under which measurements of the fluctuations of heat dis- 
sipation can be used to investigate higher order cumulants of the charge counting statistics of a 
mesoscopic conductor. 
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During the second decade of mesoscopic physics there 
has been an increasing interest in thermal phenomena. 
Several experiments succeeded in measuring properties 
of heat transport in mesoscopic samples: Experimen- 
tal tools based on the Coulomb blockade to measure ac- 
curately local teniperatures in mesoscopic samples have 
been established HQ- It has been demonstrated that 
the heat conductance through one-dimensional phonon 
modes of a micro bridge is quantized at low tempera- 
tures 0. The universal heat conductance per mode is 
given by 7r^fc^T/3/i. The same universal heat quantum is 
also found for carriers other than bosons (see |j] and Refs. 
therein). Thermopower was experimentally investigated 
for a quantum dot in the Coulomb blockade regime 
and for multiwalled carbon nanotubes 0]- The Andreev 
reflection process is ineffective for heat transport and has 
been employed to measure thermopower in Andreev in- 
terferometers. The thermopower depends on the mag- 
netic flux and shows geometry dependent time inversion 
properties j^]. 

All these experiments have in common that they inves- 
tigate mean properties averaged over time. In electrical 
transport, however, noise measurements have become a 
very useful tool to study properties of non-equilibrium 
systems which are hidden in measurements of the mean 
current (for a review see 0)- A recent milestone was the 
first successful experimental investigation of higher order 
current correlators by Reulet et al. It revealed an 

unexpected temperature dependence that was explained 
by accounting for an external measurement circuit [lTI| . 
The calculation of average currents and current fluctu- 
ations can be unified within the concept of full count- 
ing statistics that has been introduced into mesoscopic 
physics a decade ago. It was shown by Levitov and Leso- 
vik that coherent charge transfer through a two-terminal 
conductor can be seen as probabilistic process defined 
through a set of transmission probabilities There- 
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after different works investigated the statistics of other 
measurable quantities in mesoscopic structures such as 
voltage momentum 0|, and charge inside a meso- 
scopic volume [Tsj . 

In this publication we address the statistics of fluctu- 
ations of energy exchange. We investigate the statistics 
of energy flow into a subvolume of a quantum system. In 
the first part of the paper, we write the generating func- 
tion of this statistics in a Keldysh-representation [l^ . 
The "counting field" that generates moments of trans- 
ferred energy enters as a gauge field that shifts the time 
evolution of the studied subvolume relative to that of the 
rest of the system. 

In section m we apply our approach to a simple exam- 
ple which nicely demonstrates the basic features of energy 
transfer: We consider the energy dissipation from a bi- 
ased mesoscopic conductor into one of its reservoirs. The 
problem has been addressed before for a conductor at 
zero temperature [13. In [3 it was shown that dissipa- 
tion happens symmetrically, that is each reservoir dissi- 
pates half of the energy. We find that this is true only on 
average and we determine fluctuations around this rule. 
In contrast to the statistics of charge transfer through 
such a conductor, the transfer of heat is not quantized. 
In section imi we derive general relations for the statistics 
of energy exchange between a quantum system and a lin- 
ear environment. Moments of the exchanged energy are 
expressed in terms of correlators of the variable that cou- 
ples the quantum system to its environment. We apply 
these results to a mesoscopic conductor which is embed- 
ded in a macroscopic measurement circuit in section Hvl 
The energy transfer statistics is related to correlators of 
the current through the mesoscopic conductor. One ex- 
pects that the second moment of energy dissipated in a 
series resistor depends on the fourth moment of current 
fluctuations in the mesoscopic conductor. This suggests 
to use a measurement of heat fluctuations as a probe 
to study higher order current correlations in mesoscopic 
conductors. We quantify this expectation and derive the 
conditions under which higher order current correlators 
can be extracted from an energy measurement. 
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FIG. 1: Quantum system divided into two subvolumes V and 
V that are coupled via the variables jv and X^- 



Z arises naturally as one models a linear detector that 
measures the energy flux Hv into V O. The detector 
read-off r is then a linear functional of Hv , 



(r(r)) 



dts{T - t)Hv{t) 



(6) 



where the response function s{t) is causal, s{t) — for 
t < 0, and it depends on the internal dynamics of the 
detector. Moments of the detector read-off are then gen- 
erated by Z according to 



I. APPROACH 

We study the statistics of the exchange of energy be- 
tween two subvolumes V and of a quantum system. 
We assume that no degrees of freedom leave or enter V 
and that V is coupled to V only through one degree of 
freedom jv- jv shall be coupled to the variable Xv in 
V , see Fig. ^ A generalization of our approach to more 
than one coupling variable is straightforward. Our model 
Hamiltonian reads 



H = Hv + Hy+jvX^ 



(1) 



with [jV,Xy] = [jV.Hy] = [Hv,Xy] = [Hv . Hy] = 0. 

To calculate the statistics of the flow of energy into V 
in a period of time [0, r] we assume it to be in an eigen- 
state \Ey) of Hy at time t = and calculate moments 
of Hy at time t = t, defining the generating function 
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(We have set h — 1.) The average is taken over the 
initial state \Ev). The exponential exp{—iS,Hv} gener- 
ates moments of the energy in V at time r. Through 
the exponentials exp{i^Hv /2}\Ey) — exp{iS^Ey /2}\Ey) 
and {Ev\exp{i£^Hv /2} the initial energy Ey in V is 
subtracted from that such that Z generates moments of 
AHy = e'^^Hye-'"'^ - Ey, the flux of energy into V 
during [0, r], 



{{AHvr) = z 



dP - 
P—Z(f) 



(3) 



A generalization of Z that generates correlators of the 
energy flux at flnite frequency is the functional 



(4) 



The symbols T [T) denote (inverse) time ordering and 
we have collected the two source functions into a vector 
in a "Keldysh space", ^ = (C^iC )■ -2 correlates the 
energy flux at different times, 
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q=l 



(5) 



rP(r)) ( -) ^Z 



Cs(r-t) 
-(sir - t) 



(7) 



Evidently, Eq. Q reduces to (0) for the measurement of 
the time integrated energy flow during [0,t], if s{t) = 
9r{t) = e{t) - 9{t - t) [with the step function e{t) = 1 
for t > and e{t) = for i < 0]. 

In the interaction picture with respect to 
the uncoupled problem Hq ~ Hv + Hy, 
jv{t) — exp{iHvt}jv exp{—iHvt} and Xy{t) = 
exTp{iHyt}Xy exp{—iHyt}, Z takes the form 



Te 



-i f dt [Xv{t)jv(t)-i+(t)Hv] 



(8) 



We rewrite Eq. © by breaking the time-ordered prod- 
uct up into a large number of time-development expo- 
nentials for infinitesimal time steps e. Applying to each 
one of them the identity 



g-i//v«(t)g-J£iv(t)^(>(t)g*-f^v«(t) 



we find that 



(9) 



= (Te""-^ dtjv[t-rit)]Xf{t)i^^-if dtjv[t-i+(t)]x^(t)'^ 

(10) 

The generating functional of the statistics of energy flow 
into a volume V takes the form of a partition functional 
of the entire system with shifted time arguments of all 
variables in V. An analogous structure has been found 
before for the statistics of transfer of other globally con- 
served quantities, like charge 12] and momentum '14!). 
Also there the generating functional has the structure of 
Eq. 1)11) |l . The source term locally shifts the variable con- 
jugated to the measured quantity (the phase for a mea- 
surement of charge, position for a momentum measure- 
ment, and time for the energy measurement considered 
here) . 

In a path integral formulation of Eq. Q we mark flelds 
that evolve the system forward and backward in time 
with superscripts -I- and — respectively and we collect 
them in vectors like the source 1^. We may change inte- 
gration variables as jy{t — (^{t)) jyit) in the action 
corresponding to Hy if |^(i)| < 1 at all times (only then 
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Bolometer 



FIG. 2: Geometry considered in this section: A bolometer 
is attached to one reservoir of a biased mesoscopic conduc- 
tor. We calculate the statistics of heat dissipated into the 
bolometer. 



this map is bijective). The path integral takes then a 
form very similar to Eq. (|10|l . 



VjvVXy e 



xe 



-i J dtJvit-i{t))T3Xy{t) 



(11) 



We have introduced actions Sy and Sy for the uncoupled 
volumes V and V that are obtained by integrating the 
actions corresponding to Hy and Hy over all variables 
except jv and Xy respectively. T3 is the third Pauli ma- 
trix and we use the notation [jv{t—£,{t))]" = 
(a = {+, — }). The assumption \^{t)\ < 1 has the follow- 
ing physical interpretation. From Eq. JSJ and its finite 
frequency generalization it is seen that one can infer a 
probability distribution of energy flow by Fourier trans- 
forming the generating function Z{^). To resolve energy 
differences of the order of AE in this distribution one 
needs to know Z for values ^a_e — 1/AE. The condition 
\^{t)\ < 1 means then that the energy detector smears 
out the measurement over a time At > ^ab — 1/AE, 
that is, it does not attempt to measure the energy faster 
than it is allowed to by the uncertainty relation. 



II. HEAT DISSIPATION STATISTICS IN 
TUNNEL JUNCTIONS 

We first illustrate our method with a particularly sim- 
ple example: we calculate the statistics of heat dissi- 
pated from a voltage biased mesoscopic conductor into 
a bolometer attached to its right contact (see Fig. [21 . To 
this end, we express the action of a simple connector in 
circuit theory by the Keldysh Green's functions Gl,r of 
the adjacent reservoirs [l9ll20| 



5*0011 = ^ X! 



l + ir„({Gi,Gi}-2) 



(12) 



where the r„ denote energy independent transmission 
probabilities and the brackets {, } anticommutation. The 
trace includes integration over time. In the following, we 



address the zero frequency limit only. In this case, the 
time integration in Eq. (|12|) is conveniently rewritten as 
energy integration. As formulated in Eq. we obtain 
the statistics of heat dissipation into the right contact by 
shifting all observables of the right contact in time by ^+ 
on the upper Keldysh contour and by ^_ on the lower 
contour. The diagonal elements of Gr in the Keldysh 
space are left unchanged by this transformation which 
can be cast in the following form (we introduce the dif- 
ference ^ = ^+ ^ S,-)- 



1,2(1-/,) 2/,-i;- 



(13) 



The time shift introduces a rotation of the Green's func- 
tion Gfi analogous to charge counting statistics [l^ . 
fj — fj{e) denotes the energy dependent Fermi func- 
tion of contact j — L, R that depends on temperature Tj 
and electrochemical potential Hj . Evaluating Eq. p2l) in 
the zero freq uency limit we find the cumulant generating 
function [24| 

lnZ[e] = 

^ / E„ In {1 + r„/L(l - /fl)(e-^«(-^«) - 1) 
+r„/fl(l-/L)(e'«(^-^«)-l)}. 

(14) 

Cumulants of dissipated energy can now be calculated 
easily by taking derivatives with respect to ^. The mean 
dissipation without applied voltage {fj,j = 0) is for in- 
stance given by 



{AHy) = 



.dlnZ 



(15) 

For almost equal temperatures — > Tr we recover the 
mesoscopic version of the Wiedemann- Franz law for ther- 
mal conductance '2l','25|. 

We note that Eq. H14|l is similar to the statistics of 
charge transfer 0: Charge statistics is recovered by sub- 
stituting S.{£—^J'R) I— > where e is the elementary charge 
and the field x generates cumulants of charge transfer. As 
for charge statistics, there is a classical probabilistic in- 
terpretation of Eq. (|14|l : An electron of total energy e in 
channel n passes the mesoscopic conductor with probabil- 
ity r„. It then dissipates the energy e — 11 in the right 
contact. It is the energy integration in Eq. (|14|l which 
makes the statistics of heat transfer interesting. Unlike 
in the binomial charge statistics the exponent z^(e — ^n) 
assigns energy weights to each electron. Therefore, en- 
ergy transfer is not quantized. 

We now turn to a specific example: the tunnel junction 
with r„ ^ 1. In this case, the logarithm of Eq. H14|l 
can be expanded and the energy integration involves only 
elementary integrals (for equal temperatures = Tr = 
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FIG. 3: Comparison of energy dissipation AHv = Hv ~ (Hv) 
(upper panel) and charge AQ — Q — {Q} statistics (lower 
panel) for a tunnel contact in equilibrium and non-equilibrium 
due to an applied bias or temperature difference. The energy 
is normalized to temperature kTL or bias fi = hl — ^ir re- 
spectively. The product of = jir /2-k ,kTLT /2Ti and dimen- 
sionless conductance G is assumed to be large. 

T on both sides of the junction). Defining /i = /i^ — we 
find the following result for the characteristic function: 

InZK] = Gr coth ^ (kT^,0l^ - ^^) 

(16) 

sinh(7r^fcT)/7r ' 

We introduced the conductance of the junction G ~ 
^r„/27r. For comparison, we also give the statistics 
of charge transfer [23 

ln2:[x] = Gr^ coth (cos(xe) - 1) - iG/i sin(xe). 

(17) 

We observe several differences: Whereas Z[x] is strictly 
periodic, the characteristic function of heat dissipation 
Z[^] shows damped oscillations. The lacking periodicity 
is due to the unquantized transfer of energy. For charge 
transfer, the real part of Z['x\ is even in the applied volt- 
age /i in contrast to an odd imaginary part. Odd cumu- 
lants therefore change sign under voltage inversion. This 
is not the case for energy, does not depend on the 
sign of the applied voltage. Heat dissipation takes place 
in both reservoirs symmetrically regardless of the current 
direction! Odd cumulants of charge transfer (the imagi- 
nary part of Eq. 1)17(1 ) do not depend on temperature j23|. 
In contrast, we find that the asymmetry of the distribu- 
tion of dissipated heat does depend on temperature (see 
the imaginary part of Eq. H16|l ). 

Another interesting case is the statistics in the presence 
of a temperature gradient only. Analytical results are 




AHy / NG 

FIG. 4: Probability Distribution of dissipated energy into one 
reservoir of a biased mesoscopic conductor (/i ^ kT). The 
distribution is broader for tunneling junctions (F = 0.01) 
than for open point contacts (F = 0.80). It is clearly seen 
that the third cumulant of the distribution changes sign. The 
inset shows the third cumulant {{AHv)'^) as a function of 
the transparency F. (G is the dimensionless conductance and 
— fiT/2n, we assume NG to be large) 

available for Tl->Tr. We find 

^Grkn { sshfcy + ln(2) - PHkm)} = 

Grkn {-i^km - ^{km? + ^'^ikm? + ■■} 

(18) 

for energy transport and 

In Z[x] = 2GTkTL ln(2) (cos(xe) - 1) (19) 

for charge transport (we introduced the /3-function 
Piiz) = SJ^Lo(~l)"/(*'^ + '^))- Fig- El illustrates several 
limiting cases: On a log scale, it compares the charge 
and energy transfer statistics in equilibrium, and in non- 
equilibrium due to an applied voltage (see Eqs. (|16|l and 
(I17|l ) and due to an applied temperature gradient (see 
Eqs. (HHJ) and (Unj). 

In general, Eq. H14() has to be evaluated numerically. 
Fig. 01 shows the probability distribution of dissipated 
heat for barriers with various transparencies r„ = F in 
the high-voltage regime fj, ^ kT. Noise disappears in 
the ballistic limit F = 1. It is clearly visible that the 
third cumulant of the distribution changes its sign as a 
function of F. We find that 

((Ai?v)3} =.Ga.4(1-F)(1-2F)/4 (20) 

(see also the inset of Fig. 

In this entire section, we assumed that the energy dis- 
sipation in the right reservoir is measured by an ideal 
bolometer which does not act back on the mesoscopic 
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contact. A realistic bolometer is characterized by a fi- 
nite thermal conductance. Energy fluctuations in the 
right reservoir are then converted into temperature fluc- 
tuations which modulate the noise intensity of the meso- 
scopic contact. This backaction is similar to the backac- 
tion from a non-ideal current meter pjl|. 



V 



III. ENERGY FLOW INTO A LINEAR MEDIUM 

We use now Eq. to calculate the statistics of en- 
ergy flow from a system S into a linear medium. An 
example is the energy that S emits as electromagnetic 
radiation. The electromagnetic field is then the linear 
medium into which energy flows. This energy flow may 
for example be measured with a photo-detector. The dis- 
tribution of the number of photons emitted by a source 
current js can be expressed in terms of correlators of js 
[26|. In Ref. this relation has been established per- 
turbatively in a weak coupling of the photo-detector to 
the electromagnetic field. In Ref. it has been used 
to calculate the fluctuations in the number of photons 
emitted by a mesoscopic conductor. Eq. ((TT|l allows us 
to obtain these results non-perturbatively. As a second 
application of our method we calculate the fluctuations 
in the amount of heat that is dissipated in a macroscopic 
electrical circuit around a mesoscopic conductor. Such 
fluctuations could be measured with a bolometer. 



A. General relations 

In our analysis we divide space into three regions. The 
system S, a volume V into which the energy flow is mea- 
sured and a region E. We call E = E+V the environment 
of S. The subspaces are connected via two variables Xs 
and Xy in E, that couple to js in S and jV in V respec- 
tively (see Fig. O. With the Hamiltonians Hs, Hj^,, and 
Hv of S, E, and V we then have the total Hamiltonian 

H = Hs + H^ + Hv + Xsijs-js)+Xy{jv-jv), (21) 



where jv and js are sources that generate moments of 
Xy and Xs- The generating functional for correlators of 
the energy flow into V takes the form of Eq. with 



= / VXsVjs e 



-i{Ss\js]+Sii,[Xs,Xv]) 



with the actions Ss and 5^ corresponding to Hs and 
H^. S^, being the action of a linear system, is quadratic 
and depends only on the response functions of E and its 
temperature T. We characterize the response of E when 
disconnected from S (corresponding to the Hamiltonian 



FIG. 5: S is an arbitrary quantum system. The energy flow 
into the part V of its linear environment is considered in this 
section. The variables js, Xs and jv, Xy couple the three 
constituent volumes. 



H — Hs — Xsjs) by four functions 
dXsiiu) 



i?ys(w) 



djsi^) 
dXy{uj) 



djs{(^) 



Rsv{.<^) 

i?yy(w) 



dXsiuo) 



djv{uj) 
dXyiuj' 



djv{i^) 



(23) 



The volume V is described by the response to a source 
Xy coupling to jy in the absence of E, corresponding to 
the Hamiltonian Hv — jvXy, 



. djv{uj) 



(24) 



The environment's action is determined by the fluctua- 
tion dissipation theorem . The action for jv when V 
is isolated reads 



Siso[jv] = jv ® Miso ® jv- 



(25) 



The matrix multiplication ® extends over Keldysh as 
well as time indices, [A®B]{t,t") = J dt' A{t,t')B{t' ,t") 
and correspondingly for vectors like jvit)- In a 
steady state Miso{t,t') depends only on the time dif- 
ference t — t' and in Fourier representation M(w) — 
J d{t ~ t') e*'^(*"*')M(< - t') it reads 



M,so{uj) = ii?,,o(w)r3+iImi?,,„(w)Go(a;), (26) 

r (.A - ( ^ 



with the Bose-Einstein distribution N{ijj) = 
(exploj/fc^T} — 1)~^. (We define the response functions 
such that they have negative imaginary part.) 

We analyze first the particularly simple case that the 
energy flow into the entire environment to S is measured, 
V = E. Then we have Xy = js, Xs = jv and Sy = Ss- 
We rewrite Eq. by introducing a coupling matrix 



6{t-t' + (-it')) 



(28) 
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The Gaussian integrals over jy in are easily done, 
resulting in 



-i{SslIs]+SE lis ] +^|-[Js ] ) 



(30) 



with an action Se that describes the influence of the 
environment on S, 



and a source term that vanishes at ^ = 0, 

Ai^s] ^ Js(^A^-(g,js, A^-=-lT3a^3, (32) 

(33) 



In the general case V E, the environment's action is 
determined by Rss, 

Se[js] = Js®Gss(^Js, (34) 

(35) 



(31) The source term then takes the form 
I 



A, 



G 



sv ®[[Gvv + Gvv ®a^®Gvv)\ ^-G^y} 



vs- 



(36) 



Eq. I|32|l is recovered for Xy = js and Xs ~ jy , such and so is A^ . At zero temperature it takes a simple form 
that Ryy = 0, Rss — —Riso^ ^-iid Rys — Rsv = 1- even in the general case E ^V, 



B. Zero frequency 



For concrete results we focus on zero frequency corre- 
lators of the energy flow, choosing ^ = {l,—l)9r{t)£,/2 
[c.f. Eq. O]. We neglect transient effects, that is, we 
calculate only terms of leading order in the time r over 
which energy is accumulated. We may then work in 
the discrete Fourier space of r-periodic functions /, with 
Fourier coefficients /; = dt e'"'* f{t)/T and frequencies 
uji = 2'kI/t. In this representation the source term in Eq. 
(jSOJ reads 



a is now diagonal in frequency indices, 



e*'^'«/2 



(37) 



(38) 



A 



|i?ys(t^i)l' 



^Imi^^l(wO(e-'«"' 



1 



(39) 

Applied to a linear photo-detector this reproduces the 
distribution of photo-counts in response to a source js 
obtained perturbatively in |26|. To see this we substi- 
tute ^LOi — > ^„ in Eqs. (|Sn|l with l|57|l and to obtain 
the statistics of the number of absorbed quanta n rather 
than that of the absorbed energy. The further substitu- 
tion e~'^" — 1 — > Xn yields the generating function for 
factorial moments of n. To compare the resulting facto- 
rial moments of the photocount to the formulas obtained 
in i2ii, we write them in the time domain, 



(n(n — 1) • • • (n — p)) = 



di! dt" I dr' dr" ImR-^it', t")Rys(t', T')Rys{t", r")js (r')jl(r"^ 



(40) 



This is equivalent to Eq. (30) in |26j |. The detector sensi- absorbing detector modes per unit of frequency. The re- 
tivity / there corresponds to our Imi?~^, the density of tarded photon propagator Dret corresponds to our cross- 
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Zs 




Zv 







V 



FIG. 6: Electrical circuit analyzed in this section. Rs is the 
conductor that generates current fluctuations Sis- These fluc- 
tuations can be characterized by a measurement of fluctua- 
tions of the heat that is dissipated in the external resistor 
Rv- 

response function iRvs- The time-ordering of and 
jg along the Keldysh contour corresponds to the "apex" 
time order of the sources in . The expectation value 
in Eq. H40|) is taken with respect to the action 5*5 of the 
source system and the piece Se that describes the back 
action of the environment on S. 



IV. LINEAR ELECTRICAL CIRCUITS 

Electrical conductors couple to their electromagnetic 
environment - the circuit that they are embedded in - 
with the current 1$ that flows through them. To cal- 
culate the flow of energy into this environment we may 
therefore apply the formulas obtained in the previous sec- 
tion with js = I Si Hs being the Hamiltonian of the con- 
ductor. Also in this case energy is transferred by means 
of photons. These photons may be either detected by 
a photocounter or with a bolometer, that measures 
the thermal energy exchanged by means of photons. Ac- 
cording to Eq. H4U|I the n-th factorial moment of energy 
transfer is proportional to the 2n-th moment of fluctua- 
tions of the source js- This suggests that a measurement 
of the energy flow and its fluctuations may be a useful 
tool for measuring higher order correlators of electrical 
currents whose measurement poses an experimental chal- 
lenge. While experimental techniques for the measure- 
ment of the variance of electrical current fluctuations are 
by now well developed, so far only one experiment has 
been successful in measuring higher order current corre- 
lators. In lO] the measurement of the third moment of 
current fluctuations through a tunnel barrier has been 
reported. Eq. H40|l suggests that the fourth moment of 
current fluctuations can be inferred from the variance of 
the heat produced by them. 

To quantify this we apply now the relations obtained 
in the previous section to the electrical circuit depicted in 
Fig. El We assume all resistors in the circuit to be macro- 
scopic and linear, except i?s, that plays now the role of 
the system S. Fluctuations of the amount of heat that 
is dissipated in Ry as measured with a bolometer are by 
virtue of Eq. related to fluctuations in the electrical 



current through Rs- For that the volume V should be 
chosen the resistor Ry- The current /y flowing through 
the volume V is coupled to the variable Xv '"^^t of 

the environment E. This variable Xv the time integral 
of the voltage Vy over i?y, Xv(t) = fdi'Vyit')- The 
current Is through Rs similarly couples to the environ- 
ment variable xs(t) = J* dt' Vs{t')- Again we introduce 
external sources ly and Xv- The circuit Fig. |H|is then 
mapped onto the general model of Fig. |5l with the choice 

=_ Xy, = Xy. = Xs, jv = W, jv = 
js — Is and js = Is - We need the response functions of 
the circuit Fig. without the conductor Rs, 



Rss 



dx& 



dis 



Iv 



1 

iuj 



1 



Zy + Ry 



Rys — Rsv — 



dxs 



dh 



1 



ZR\ 



Is 



R 



vv 



9x\ 



dh 



Is 



1 

iuj 



iuj Z -\- Z\ 



1 



Ry 



Zv 



Z 



dh 



dx\ 



IUJ 

isolated Ry 



, (41) 



(42) 



(43) 



(all quantities here are frequency-dependent). Note, that 
any four-terminal circuit connecting the conductor Rs 
with Ry can be modeled by the three resistors Z , Zs, and 
Zy - The resulting zero-frequency energy flow statistics 
takes a particularly simple form in the limit of an infinite 
Z . Then the current flowing through Rs is directly fed 
into Ry and the source term in the generating functional 
is given by Eq. with 



A 



Re Ry {uji) 



This is a generalization of Eq. (|39|1 to finite temperature, 
allowing for emission of energy quanta from V as well as 
absorption. The transferred quanta have predominantly 
energies that are smaller than the inverse of the RC-time 
'Trc of Ry. For finite Z thermal current fiuctuations 
{5P)th ~ kT/Z created in Z mix into the fluctuations of 
Is- As a consequence Eq. (|44|l is then only valid in the 
regime {TRckTf < \Z/Ry\- 

The formulas for the statistics when the mixing in of 
thermal fluctuations occurs are more complicated. To 
make further progress we assume that the frequency dis- 
persion of the measured current correlators is negligible 
on the scale 1/trc, ((Ilg^i -^s('^g))) ~ 27r(5(X;^=i Wg)Cp 
for LOq ^ I/trc (here ((• • • )) denotes irreducible, or cu- 
mulant correlators). For a mesoscopic conductor this is 
satisfied if I/thc is smaller than the voltage applied to 
the conductor. In this limit we find that 
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{AHy) = 2r ^R,ff{u;)Q 
{{{AHyf)) = At- 



zvr 



-1 2 










io 2^ - 



(45) 



(46) 



with Reff — Rei?y|Z/(Z + + Rv)\^- We have assumed that i?e//(0) = and Reff{uj) — > for a; > I/trc such 
that the DC-component (that wiU be avoided in experiments) and the contributions from frequencies oj > 1/tjic 
in Eqs. (|45l) and H4()l) are neghgible . We conclude that the variance of the heat dissipated in Ry depends on the 
fourth cumulant of current fluctuations in the mesoscopic conductor, as one would expect. There are, however, 
also contributions from lower order current correlators. The environment circuit has a finite response time tuc ^^nd 
therefore effectively averages fluctuations over that time. As a result, the energy fluctuations become dominated by 
the lowest order current correlator in the limit of a long trc, when higher order current correlators become negligible 
and the statistics becomes Gaussian. In this limit of "narrow-band detection" the statistics of energy transfer is 
negative binomial '27|. Only deviations from this encode non-Gaussian current correlations. In order to see them 
the statistics of charge flow through the conductor during tjic has to be strongly non-Gaussian. This is the case if 
trcI /s, the mean number of transmitted electrons in that period, is small. 

More concretely, for a measurement of the fourth cumulant C4 of current fluctuations in Rs one would want the 
first term in Eq. (|4(j|) to be dominant. An estimate of the three summands in Eq. H46|l for a tunnel contact with mean 
current / (assuming that Reff {to) ~ Reff for uj < I/thc), 



{{{AHyY)) « —Rl Ggel 1 



' RC 



TRC I _^ T|.pfcrmin{fcr, 1/trc} 

e GqReff 



(47) 



(with the conductance quantum Gq — /2ti) confirms 
this. The variance of fluctuations of the heat flux into 
Ry is a direct measure for the fourth cumulant of cur- 
rent fluctuations in a tunnel contact if t^c^ je <^ 1 and 
TRckT -C GqReff. The origin of the first condition has 
been explained qualitatively above. The second condi- 
tion ensures that the fourth cumulant is visible on the 
thermal background. Back action effects of the measur- 
ing resistor on the measured tunnel contact (conductance 
G) are avoided if additionally RssG <C 1 These re- 
quirements are rather restrictive, but can in principle be 
met in experiments. For a practical implementation, the 
first condition tjicI / e <^ 1 may be relaxed to trcI — 1 
by measuring the voltage dependence of {{{AHy)"^)) and 
extracting the term that is linear in /. 

V. CONCLUSION 

We have presented a theory for the statistics of en- 
ergy exchange between coupled quantum systems. As 
an application we have calculated the statistics of energy 
dissipation into the leads connected to a mesoscopic con- 
ductor. General and exact expressions can be obtained 
for the energy flow from a quantum system into a linear 



environment. We have used these results to calculate the 
moments of the heat dissipated in a linear circuit around 
a mesoscopic conductor. They have been expressed in 
terms of moments of the current fluctuations produced 
by the conductor. As one expects, the variance of this 
dissipated heat depends on the fourth cumulative mo- 
ment of current fluctuations produced by the conductor. 
Heat fluctuations may therefore serve as a tool to detect 
the fourth cumulant, whose direct measurement is dif- 
ficult. We have analyzed the conditions under which a 
measurement of heat fluctuations can reveal higher order 
cumulants of the charge counting statistics. 
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